DESY 99-037 
August 1999 
|hep-ph/9908289 



PREDICTIONS OF ZFITTER v.6 FOR 
FERMION-PAIR PRODUCTION WITH ACOLLINEARITY CUT 



P. CHRISTOVA^^ 
Faculty of Physics, Bishop Konstantin Preslavsky Univ., 
Shoumen, Bulgaria 
and 

Lab. for Nuclear Problems, Joint Institute for Nuclear Research, 

Dubna, Russia 
E-mail: penchris@nusun.jinr.ru 

M. JACK, S. RIEMANN, T. RIEMANN 
DESY Zeuthen, Platanenallee 6, 
D- 15738 Zeuthen, Germany 
E-mails: jack@ifh.de, riemanns@ifh.de, riemann@ifh.de 



ZFITTER is a Fortran package for the description of fermion-pair pro- 
duction in e^e~-annihilation. We report on results of a rederivation of 
the complete set of analytical 0{a) formulae for the treatment of pho- 
tonic corrections to the total cross-section and the integrated forward- 
backward asymmetry with combined cuts on acoUinearity angle, ac- 
ceptance angle, and minimal energy of the fermions. Numerically, the 
following changes result in ZFITTER v. 6.11 compared to ZFITTER v. 
5.20/21: (i) at the Z resonance - numerical changes are negligible; (ii) 
at LEP 1 energies off-resonance - corrections amount to at most few 
per mil; (iii) at LEP 2 energies - corrections amount to one per cent 
or less. Thus, the predictions for LEP/SLC data remain unchanged 
within the actual errors. 



'^Based on talks presented at ECFA/DESY Linear Collider Project meetings held at Frascati, Nov 
1998, and Oxford, March 1999, at the LEP 2 Mini- Workshop held at CERN, March 1999, and at 
the Workshop for a Worldwide Study on Physics and Experiments with Future Linear e'^e" Colliders 
held at Sitges/Barcelona, April 1999. 
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1 Introduction 



The Fortran package ZFITTER |],|,|@J|J|] is regularly used for the analysis of LEP 
data since 1989 for the reaction: 

e+(fc2) + e'ih) ^ f{p2) + fipi) + {n^){p). (1) 

Since the systematic description of ZFITTER v.4.5 (19 April 1992) [0 a series of improved 
versions has been released and a nearly complete collection of them may be found in 
the world wide web Pfl. Most important updates of recent years are related to the 
treatment of higher-order corrections. In 1995, this was documented systematically in 



I0| . Later improvements are described in [jl|. The ZFITTER versions 5.20/21 |Tl| were 
quite recently followed by versions v.G.nm, beginning with v. 6. 04/06 [Q. 

Traditionally we aimed at an accuracy of ZFITTER at LEP 1 energies of the order 
of 0.5 %. The successful running of LEP 1 together with the precise knowledge of the 



beam energy however makes an even higher precision necessary [13]: We expect for 
the final measurements relative errors for total cross-sections and absolute errors for 
asymmetries of up to 0.15% at the Z peak and of up to 0.5% at y/s = Mz ± several 
GeV; aiming from the theoretical side ideally at a tenth of these values for the errors 
of single corrections, we estimate limits of 0.015% and 0.05%, respectively |. Similar 
claims may be found in [|T^ . 

First applications of ZFITTER at energies above the Z resonance have become rel- 
evant since data from LEP 1.5 and LEP 2 are being analysed. There is also rising 
interest in applications for the study of the physics potential of a Linear Collider oper- 
ating at a/s > 500 GeV |]T6[. At the higher energies, again one generally expects final 



experimental accuracies of up to 0.8% [T^. So, one should aim at theoretical accuracies 
of single corrections of up to 0.1 %. The excellent precision of ZFITTER at the Z peak, 
however, does not automatically guarantee a sufficient accuracy at higher energies, es- 
pecially since the hard photonic contributions, including higher-order corrections, are 
no longer suppressed; see Figure |I] and compare to Figure |. 

In 1992, a comparison of ALIBABA v.l (1991) and ZFITTER v.4.5 (1992) showed 
deviations between the predictions of the two programs of several per cent |jl8[; we 
show one of the plots of that study in Figure These deviations were observed only 
above the Z peak and only when an acollinearity cut on the fermions was applied; the 
agreement was much better without this cut. We repeated the comparison in 1998 
with ALIBABA v.2 (1991), TOPAZO v.4.3 (1998) 0,0, and ZFITTER v.5.14 (1998). The 



outcome was basically unchanged compared to 1992 as may be seen in Figure 3 of 
reference \T3\. 



It is well-known that deviations of up to several per cent may result from different 
treatments of radiative corrections. Recent studies |^ claim for the case of Bhabha 



scattering that an accuracy of 0.3% for 0{a) corrections and of 1% for the complete 
corrections has been reached at LEP 2 energies as long as the radiative return to the 
Z peak is prevented by cuts. Similar conclusions were drawn in for fermion-pair 
production. The radiative return is prevented if > Mz, i.e. if R = s' / s > M^/ s. 

'^This might not even be sufficient if the Giga-Z option of the TESLA project will be realised (see e.g. 
jOj), with a factor of 10 or 100 more Z bosons produced than at LEP 1. 
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Figure 1: Cross-sections and forward-backward asymmetries for muon-pair production with different 
acollinearity cuts. 



Our figures contain predictions with an acollinearity cut of ^acoi < 10°, 25°. Tliis 
corresponds to an s'-cut of roughly R > -Rgacoi — 0.84,0.64 and y/s > 100 GeV, 114 
GeV, respectively (see also Table |A.l] in the Appendix). The resulting suppression of 
the radiative return at these and higher energies can nicely be observed in Figure |I|. 
But even if the radiative return is prevented, the influence of hard photonic corrections 
will be much larger at higher energies than it is near the Z resonance where hard 
bremsstrahlung is nearly completely suppressed. Figure ^ demonstrates that different 
portions of hard photon emission lead to nearly identical cross-sections unless the region 
is reached where even soft photon emission is touched (lowest lying curve). 

In this paper, we report on a recalculation of the photonic corrections with acolli- 
nearity cut in the ZFITTER approach. In Section 0, we explain changes of flags, sub- 
routines, kinematics, and numerics of ZFITTER. The basic differences to the simpler 



s'-cut are explained in subsection |2.2| . Section y contains few comparisons with other 
programs and Section ^ a summary. 

2 Photonic Corrections with AcoUinearity Cuts 

Perhaps we should remark here that the photonic 0{a) corrections in ZFITTER remained 
basically untouched since about 1989. For the s'-cut, ZFITTER relies on duplicated 
analytical calculations p3| , p4| ,p|,^ , and numerical comparisons showed the reliability of 



the predictions at LEP energies; see e.g. [O] for LEP 1 and [EH] for LEP 2. Concerning 



the acollinearity cut, the situation is different. The corresponding part of ZFITTER |26] 



was never checked independently, and is not documented]^. So, we started a complete 

A collection of some formulae related to the initial-state corrections (and its combined exponentiation 
with final-state radiation) for the angular distribution may be found in . 
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Figure 2: Muon-pair production cross-sections at LEP 1 energies with different acollinearity cuts. 
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Figure 3: IVluon-pair production cross-section ratios ALIBABA v.l (1990) versus ZFITTER v.4.5 (1992). 
An acollinearity cut was applied with = 25°. 
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revision of the 0{a) contributions with acolhnearity cut. First results were pubhshed 
in |TH] where we reported that some deviations from ZFITTER were found for initial- 
state radiation. We also reported that these deviations could not explain the differences 
observed when comparing with ALIBABA, while the treatment of higher order corrections 
evidently was of much higher influence in this respect. Slightly later it was observed 
in |15| that the perfect agreement of many predictions of ZFITTER v. 5. 20 and TOPAZO 
v.4.3 at LEP 1 energies of about typically 0.01% could not be reproduced when an 
acolhnearity cut was applied and the initial-final state interference was taken into 
account. This second puzzle could be resolved by our recalculation; see section [Z^ . 

By now we have a complete collection of the analytical formulae for the 0{a) correc- 
tions. The corresponding Fortran package is acol . f . We merged package acol . f with 
photonic corrections for the integrated total cross-section and the integrated forward- 
backward asymmetry (with and without acceptance cut) into ZFITTER v.5.21, thus 
creating ZFITTER v.6. 04/06 |T^ onwards. The angular distribution will be available 
in v.6. 2 onwards. The remarkably compact expressions for the case that no angular 
acceptance cut is applied are published 0. A complete collection of the analytical 
expressions is in preparation [^ . 

2.1 Flags and subroutines in ZFITTER 

In ZFITTER the calculation of cross sections and asymmetries is dealt with by subroutine 
ZCUT which either calls subroutines SCUT (angular acceptance cut applied) or SFAST 
(shorter formulae; no angular acceptance cut). In the latest releases, ZFITTER v. 5. 20 
up to ZFITTER v.6. 11, the following flag settings in subroutine ZUCUTS are of relevance 
here: 

• ICUT = -1: cut on the invariant mass of the fermion pair, s' 

• ICUT = 0: cuts on acolhnearity and minimal energy of the fermions and on the 
acceptance angle of one fermion 

• ICUT = 1: cuts on s' and on the acceptance angle d of one fermion 
Since ZFITTER v.6. 04, we introduced two new values of this flag: 

• ICUT = 2: new coding of case corresponding to ICUT = 0, but without acceptance 
cut [01; 



ICUT = 3: new coding of case corresponding to ICUT = |[28| . 



In order to maintain compatibility with earlier releases, the branch ICUT=0 is retained; 
it calculates with the old coding of the final-state corrections (flag IFUNFIN = 0), while 
all other branches use the newly corrected final-state contributions (flag IFUNFIN = 1). 

All the corrections related to the acolhnearity cut are called from essentially only six 
subroutines/functions of ZFITTER; see Figure ||. The subroutines SHARD and AHARD, or 
respectively FCROS and FASYM, call the initial-state and initial-final state interference 
corrections, while FUNFIN calculates final-state corrections (for the case of common 
exponentiation of initial- and final-state soft photon corrections). 

For a detailed description of the Fortran package ZFITTER we refer to 'M. 
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Figure 4: Logical structure of subroutine ZCUT. 



2.2 Photonic corrections with acolhnearity cut: Kinematics 

The phase-space parameterisation 



The phase-space parameterisation derived in |^ is used. A three-fold analytical in- 
tegration of the squared matrix elements had to be performed. The last integration, 
that over R = s'/s, is then performed numerically. The Dalitz plot given in Figure |^ 
may help to understand the relation between a kinematically simple s'-cut and a more 
involved acolhnearity cut. The variable shown besides R is x, the invariant mass of 
(/ + 7) the cms. As Figure ^ shows, we have to determine the cross-sections in three 
phase-space regions with different boundary values of x at given R: 



da 



III. 



dcos^ 

Region I corresponds to the simple s'-cut. 
The integration over R extends from R 



dR dx 



to 1, 



da 



dRdxd cos § 



(2) 



R„ 



Rf 



sin2(C-/2) 



(3) 



/ 1- Re cos^ (6^1^/2) 

The soft-photon corner of the phase space is at i? = 1. Thus, the additional con- 
tributions related to the acolhnearity cut are exclusively due to hard photons. The 
boundaries for the integration over x are, for a given value of R: 

1 



X 



max,min 



[l-R) (1±A), 



(4) 
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where the parameter A = A{R) depends in every region on only one of the cuts apphed: 



with: 



An 
Am 



'1 



l + R- 2Re 



1, 



l-R 



\ 



Rf 



'min 



i 

^ 1 - sin(grof/2) 



2E, 

Re = 



1 + sin(0S/2) 



(5) 
(6) 

(7) 



(8) 
(9) 

(10) 



Here, rrif and Emm are the final-state fermions' mass and a cut on their individual 
energies in the cms. 




Figure 5: Phase-space with cuts on acollinearity and energy of the fermions. 
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Mass singularities from initial-state radiation 

In the last section we saw that the Dalitz plot in Figure ^ is independent of the 
scattering angle cos Here we will show that the integrations in regions II and III are 
nevertheless crucially influenced by cost?. 

The squared matrix elements contain, from initial-state radiation, the electron 
(positron) propagator, and these terms are proportional to flrst and second powers 
of 



1 1 1 



Zi(2) ih(2) - V? - f^l 2k^2)P 
1 



Ai(2) ± -Bcosy?^' 



with 



A(2) = |(1 - R){l±Pocos^cos9^), (12) 
B = ^{l- R)Posin^sme^, (13) 



and 



/?o = Vl-4m2/s, (14) 

cos6'^ = — — ^]^5^ 
2 J A2Ap 



Ai = (l-x)s, (16) 
\2 = {x + R)s, (17) 



Ap = (1 - R)s. (18) 

For the calculation of the initial-state corrections we may neglect flnal-state mass ef- 
fects. 

The flrst analytical integration is over (p^, the photon production angle in the (/+7) 
rest system [p^] : 



dT j^^ da ^^g^ 

dRdxd cos "d Jo dRdxd cos 'dd^p^ 

It is important to take care of the electron mass me, e.g. in the following contribution: 



Zi{R, cos 'd,x, if .y) \fC., ' 



Ci = s^aiX^ — 2sbiX + Ci, (21) 
a, = s\z^-Rr]l), (22) 

h = s^[Rz, (1 - ^.) - ^ ^(1 - R) Vl], (23) 



Zl(2) 



1 T Po cos 1 ± Po cos {} 
1" -fL 



Ami 



(24) 
(25) 

(26) 



Be the next integration that over x with hmits (H). Consider e.g. the following integral: 



with 



/°(i?, cos^?) 



sdx 

7^ 



In 



^iC^ + isaiX — hi 



(27) 



max .mm 



s\l-Ry iy,±Az,y + Ril-A')r^', 



(28) 



ym 



I T Po cos {} 1 ± Po cos {} 

z R z 



(29) 



and 



SttiX 



R){yi±Azi) + 0{v, 



(30) 



We are interested in the limit of vanishing electron mass rrie for the subsequent inte- 
grations. In this limit there will occur zeros in arguments of logarithms like that in 
(|27| ) at four locations defined by the conditions: 



(yi ± Azi 



0, 



1,2. 



(31) 



These zeros appear as functions of cos-i? with parameters R and A = A{R) at certain 
values cos '(9 = cf (i = 1, 2): 



1-R + A{R){1 + R) 
1 + R + A{R){1- R) 
1- R- A{R){1 + R) 
1 + R-A{R){1-R)' 

-ct < 0. 



> 0, 



(32) 

(33) 

(34) 
(35) 



The relations C2 > C2 and cf < c{ are also fulfilled. In the course of the integrations 
different analytical expressions have to be used in different kinematical regions when 
neglecting rrie wherever possibleFI. As a result of all that we have to cut the remaining 



°In logarithmic contributions of the type Lg = ln(s/TOg) and ln(l ± /3ocosi9) from coUinear photon 
emission, the electron mass has to be taken into account. 
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phase space for cost? (at fixed R and for given i) into three different regions. It is 
region I (with the s'-cut) where the conditions ( P^D to ( P^D become triviah 

4(J) = -C2 (/) = cr(/) = -c+(/) = 1. (36) 

There, only one case has to be considered (—1 < cos?? < 1). In regions II and III, the 
differential cross-section da / {dRdcosd) is a double-sum over i = 1,2 and has different 
analytical expressions for each combination of the kinematical ranges defined by (|32[) 
to (§3). 

The final result for e.g. J° after integration over x setting me = becomes {i = 1, 2): 
(i) for I cos^9| < \c~\ with yi ± Azi > (case "i + +"): 



J_ln I + 

szi \yi-Azi 



(37) 



(ii) for |cj I < I cosi9| < \cf\ with yi + Azi > and yi — Azi < (case "z + 

^»^(?/» + ^^»)(^^» -l/») 

i?2(l-/?2pOs2 79) 

(iii) for I cos'd\ > \cf\ with yi ± < (case "i ' 

jO = -^In ^ ^^ + ^^* 





SZi 



^Vi - Azi^ 



(38) 



(39) 



c. 
d. 



"1 + +'' 



)) . 

) 

1 



One may show that the resulting number of cases for the angular distribution, depend- 
ing on the value of cos d with respect to cf and on R, is at most four in regions II and 
III. These are for cost? > with the abbreviations given in eq. (0) to (^): 

a. " 1 + +" combined with " 2 + +" ; 

b. "1 + -" combined with "2 + 
combined with "2 + 
combined with " 2 + -|-" . 

For cos -i? < cases c. and d. are exchanged by: 

a. " 1 + V combined with " 2 + V ; 

b. " 1 + -" combined with "2 + -'' 
c' " 1 + -" combined with "2 + +'' 
d.' " 1 + +" combined with "2 " . 

For region I (the s'-cut) only case b. is possible. This simplification follows from A = 1; 
see above. 

When integrating over cosd within acceptance cut boundaries ±c, the distinction of 
different regions in phase space has to be repeated where the cut-off c now plays the role 
of cos d. Depending on the relative position of c with respect to the values cf we have 
to integrate over different expressions of the angular distribution (indicated above by 
the distinction of to d., or a. to d.' respectively) in the remaining cosd-R phase 

space. One finally gets at most four different analytical expressions for a!}?'^'^{R; c, A) 
and six for ap°'^'^{R; c. A) in different regions of phase space. This is because of symmet- 
ric cancellations when integrating over cos^: aip^'^ = Jl,, d cos da^°'^'^/d cos while 
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(j^ard ^ {^^Q — /°^^ dcosdda^"-'^''- /dcos-d. It is the additional occurance of c = in the 
definition of a^'g'^ that leads to more cases. 

If no acceptance cut is applied, c = 1, only one case remains for a^"'''^ = a{l)—a{—l) 
in regions II and III - because then — 1 < cf < 1 (see cases d. and d' above) - while 
for a^^'^ two of them are left because the additional integrated contributions from 
cos '(9 = depend on whether > or not = — cf): cf'^b^ = cr(l) — 2cr(0) + cr(— 1). 
The conditions (|3^) and (p4D are fulfilled for cos d = Q with 

so that, depending on A(_R), one or the other analytical expression has to be used. 
This was practised in 0. 

One can check that the integrated results 0"^°^'^ are continuous when c cf while 
da^"''^'^ / dcosd can be regularised at cosd = cf taking the exact logarithmic results in 
rrie for the integrals. 

One can also reassure oneself that the contributions proportional to the Born cross- 
section cr° and Born asymmetry ApB are (anti) symmetric respectively as it should be 
for the one loop corrected initial-state results. We will come to this in the next section. 

The phase-space splitting discussed above has also an infiuence on the initial-final 
state interference corrections since there the initial-state propagators with appear 
linearly. They will not contribute in the final-state contributions so that the phase- 
space splitting is not necessary there. 



2.3 Initial-state radiation 

We systematically compared numerical predictions of ZFITTER v. 5. 20 and ZFITTER 
v.6.11 with default fiag settings. Version 5.20 was used as released, while version 6.11 
was prepared such that the changes due to the recalculation of initial-state corrections, 
final-state corrections, their interferences, and the net effect could be isolated. We 
begin with a study of the changes related to initial-state radiation. 

For (Jt, the changes are at most one unit in the fifth digit at LEP 1 energies and thus 
considered to be completely negligible. In Table |I] we show the corresponding shifts 
of predictions for Aps for two acollinearity cuts, 6'acoi < 10°, 25°, and three different 
acceptance cuts, 'dace = 0°,20°,40°. The changes are also less than the theoretical 
accuracies demanded. We checked that our numbers for ICUT = agree with the 
ZFITTER predictions shown in Tables 26 and 27 of |jl5|. 



In Figure |^ the ratio of Ut from ZFITTER v.6.11 and v.5.20 and in Figure |^ the 
difference oi Aps are shown in a wider energy range. While at energies slightly above 
the Z peak the differences of the predictions show local peaks, at LEP 2 energies and 
beyond they are negligible for ax and amount to only 0.1% - 0.2% for ApB- The 
peaking structures disappear at energies for which the radiative return is prohibited by 
the cuts. In dependence on the acollinearity cut, this happens for energies s > s™**^. 
The s™" is calculated in Appendix 0. 
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with ^acol < 10° 


^acc 




- 1.8 




+ 1.8 


+3 


0° 


-0.28462 


-0.16916 


0.00024 


0.11482 


0.16063 




-0.28453 


-0.16911 


0.00025 


0.11486 


0.16071 


20° 


-0.27521 


-0.16355 


0.00032 


0.11141 


0.15602 




-0.27506 


-0.16347 


0.00035 


0.11148 


0.15616 


40° 


-0.24230 


-0.14398 


0.00045 


0.09881 


0.13868 




-0.24207 


-0.14386 


0.00050 


0.09893 


0.13891 


with ^acol < 25° 


^acc 


M,-3 


- 1.8 


M, 


+ 1.8 


A'4 + 3 


0° 


-0.28651 


-0.17051 


-0.00043 


0.11292 


0.15680 




-0.28647 


-0.17049 


-0.00043 


0.11293 


0.15682 


20° 


-0.27727 


-0.16499 


-0.00038 


0.10942 


0.15201 




-0.27722 


-0.16497 


-0.00037 


0.10944 


0.15204 


40° 


-0.24452 


-0.14549 


-0.00027 


0.09675 


0.13449 




-0.24445 


-0.14545 


-0.00026 


0.09678 


0.13454 



Tabic 1: Comparison of ZFITTER v. 6. 11 (first row) with ZFITTER v. 5. 20 (second row) for the muonic 
forward- backward asymmetry with angular acceptance cut (6'acc = 0°,20°,40°) and acollinearity cut 
(6'acoi < 10°, 25°); = 91.1867 (GeV). The initial-final state interference is switched off and only 
initial-state radiation is corrected. 



For the case of initial-state radiation, we were able to trace back the reason of the 
numerical inaccuracies related to the acollinearity cut of ZFITTER below version 6. It is 
the result of leaving out a certain class of non- logarithmic, simple terms of order 0{a). 
For (Tr, polynomials proportional to cos-i? (and their integrals) are concerned, and for 
cTpB polynomials of the type (a -|- focos^-i?) (and their integrals). At first glance the 
corresponding contributions seem to vanish for symmetric acceptance cuts. But this is 
not the case! As we explained in Section |2]^, the cross-section formulae lose the usual 
simple symmetry /anti-symmetry behaviour under the transformation cos "(9 ^ ( — cos "i?) 
in regions II and III of the phase space since different analytical expressions may be 
needed depending on the location of the parameters cf describing the solutions of (pl|). 
Then, the symmetry behaviour as a function of cost? is "hidden" since different regions 
contribute differently to the net result. 



2.4 Initial- final state interference 

The ZFITTER v. 5 predictions of photonic corrections from the initial-final state inter- 
ference also receive modifications due to the recalculation for the versions v. 6. The 
explanation given in Section |2.3| for the case of initial-state radiation is also applicable 
for a part of the deviations here. The codings for the initial-final state interference also 
show additional deviations in the hard photonic corrections and the resulting numerical 
differences are much larger. 
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CT^ [nb] with ^acoi < 10° 


[nb] with 6'acoi < 25° 


^acc 




- 1.8 




+ 1.8 


+3 


^acc 


M,-3 


- 1.8 




+ 1.8 




Z6 0° 


0.21928 
n 91 779 

-7.16 


0.46285 
-4.41 


1.44780 
1 4477(S 

-0.03 


0.67721 
n (S7RQR 

U.U 1 Ot/O 

+2.60 


0.39360 
+3.27 


Z6 0° 


0.22328 
n 9999S 

-4.51 


0.46968 
n 4fi8'^fi 

u .^uoou 

-2.82 


1.46598 
+0.03 


0.68688 

U . UOO 1 u 

+ 1.86 


0.40031 
n 401 9R 
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Table 2: Comparison of ZFITTER v. 6. 11 (first row) with ZFITTER v. 5. 20 (second row) for muon-pair production cross-sections with angular acceptance cuts 
(^acc = 0°,20°,40°) and acollinearity cut (^acoi < 10°, 25°). First row is without initial-final state interference, second row with, third row the relative effect 
of that interference in per mil. Final-state treated as in v.5.20. 




Figure 6: Ratios of muon-pair production cross-sections predicted from ZFITTER v. 6. 11 and v. 5. 20 without 
and with acceptance cut and with three different acollinearity cuts: 6'acoi < 10°, 25°, 90°; Emin — 1 GeV; 
programs differ by initial-state radiation. 



In Table ^ we show the shifts of predictions for muon-pair production due to the 
initial- final state interference for two choices of the acollinearity cut. Table ^ shows 
the corresponding effects for the forward-backward asymmetry. The two tables are 
the analogues to Tables 37-40 of ]15|], where TOPAZO v.4.3 and ZFITTER v.5.20 were 



compared. At the Z peak, the predictions for the influence of the initial-final state 
interference from ZFITTER v.5.20 and ZFITTER v.6.11 deviate from each other only 
negligibly, with maximal deviations of up to 0.015%. At the wings, the situation is 
quite different; we observe deviations of up to several per mil for cross-sections and up 
to a per mil for asymmetries. The deviations between the two codings decrease if the 
acollinearity cut is weakened. 

In Figure § the corresponding shifts are shown for ut (relative units) for 6acc = 
0°, 40°, and in Figure | with 9acc = 0°, 40°, for ApB (absolute units) in a wide range of 
energies. For the cross-sections, the deviations may reach at most up to 1% at LEP 2 
energies, while for asymmetries they stay below 0.5% there. Both shifts are more than 
the precision we aim at for the theoretical predictions. 

2.5 Final-state corrections 

For the case of final-state radiation, common soft-photon exponentiation together with 
initial-state radiation is foreseen in ZFITTER. For an s'-cut, ZFITTER follows [^]. As 
may be seen from ||2^ (for the angular distributions) or from (for integrated observ- 



ables), the predictions for common soft-photon exponentiation include one additional 
integration, namely that over the invariant mass of the final-state fermion pair at a 
given reduction of s into s' after initial-state radiation^. This additional integration is 



*With acollinearity cut, there remains some arbitraryness in the choice of the region with exponenti- 
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Table 3: Comparison of ZFITTER v. 6. 11 (first row) with ZFITTER v. 5. 20 (second row) for the muonic forward-backward asymmetry with angular acceptance 
cuts (^acc = 0°,20°,40°) and acollinearity cut (^acoi < 10°, 25°). First row is without initial-final state interference, second row with, third row the relative 
effect of that interference in per mil. Final-state treated as in v.5.20. 
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Figure 7: Differences of muonic forward-backward asymmetries predicted from ZFITTER v. 6. 11 and v. 5. 20 
without and with acceptance cut and with three different acollinearity cuts: ^?acoi < 10°, 25°, 90°; Emin — 1 
GeV; programs differ by initial-state radiation. 



called from subroutine FUNFIN and is treated in ZFITTER in a mixed approach. For 
not too involved integrands, the integration was performed analytically, while the inte- 
gration of two logarithms is being done numerically using the Lagrange interpolating 
formula for the integrand (subroutine INTERP and functions FALl and FAL2). A lattice 
of 20 points is used with subroutine INTERP, and defining a more dense lattice gave no 
improvements. In ZFITTER v. 5, we found a wrong sign of the term D2 in variable SFIN 
in subroutine FUNFIN and a wrong definition of the argument of function FAL2. The 
former contributes to 0"^, the latter to (Tfb- Additionally, we checked the numerical 
stability related to the numerical integration. We see also no problem related to a 
neglect of some cos-i? dependent terms for ax in subroutine FUNFIN. 
The variables mentioned above are defined in 



For LEP 1, the numerical outcome of our minor improvements is shown in Table 
I (for 6'acoi < 10°, 25°) for ApB at several different acceptance cuts: 'dace = 0°,20°,40°. 



Again, our numbers for ICUT = agree with those shown in Tables 26 and 27 of JTS 
All the changes are though visible, but negligible. For the cross-sections, the differences 
are completely negligible and not tabulated here. 

In Figure |10|, the corresponding shifts from ZFITTER v. 5. 20 to v. 6. 11 are shown for 
(Tt (relative units) and for ApB (absolute units) in a wide energy range for the case 
without acceptance cut. We see that the deviations are also negligible in the wide 



ation. We did not change what was realised in ZFITTER v. 5: Initial-state radiation is exponentiated 
for R > max(i?£;, i?^) and the final-state radiation, at given R, for R' > laax^Rmm/ R, Re)- A 
preferred condition might be R' > inax{Rmin/R,RE,Ri)- In this case, for Rmin/R < R^, the non- 
exponentiated hard photonic corrections from region III would have to be left out in order to avoid 
double counting. 
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Figure 8: Ratios of muon-pair production cross-sections predicted from ZFITTER v. 6. 11 and v. 5. 20 without 
and with acceptance cut and with three different acollinearity cuts: 6'acoi < 10°, 25°, 90°; Emin — 1 GeV; 
programs differ by initial-final state interference. 



energy range, never exceeding 0.01% for the cross-section and 0.1% for the asymmetry. 
If an acceptance cut is apphed, the changes are yet smaller. 



2.6 Net corrections 

Finally, we want to show the resulting effects of the photonic corrections discussed 
in the foregoing sections. We have to distinguish two different approaches to data. 
Sometimes experimentalists subtract the intitial-final state interference contributions 
from measured data, and sometimes the interference effects remain in the data sample. 

The net corrections without initial-final interferences are negligible for the cross- 
section. At LEP 1, they are shown for the muonic forward-backward asymmetry in 
Tables |. 

For a wider energy range, they are shown in Figures 0-IT^. Again, at LEP 2 



energies, the changes are below what is expected to be relevant. 

We saw that the corrections to the numerical output from ZFITTER with acollinear- 
ity cut increased when the corrected initial- final state interference is taken into account. 
The resulting net corrections for the muon production cross-section and the forward- 
backward asymmetry at LEP 1 are shown in Table ^ and in a wider energy range in 
Figures to |1^. The numerical effects are dominated by the initial-final state inter- 
ference and never exceed 1% at LEP 2 energies. At LEP 1 they are much smaller; see 



the discussion in Section 12^4 . 
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Figure 9: Differences of muonic forward-backward asymmetries predicted from ZFITTER v. 6. 11 and v. 5. 20 
without and with acceptance cut and with three different acollinearity cuts: ^?acoi < 10°, 25°, 90°; Emin — 1 
GeV; programs differ by initial-final state interference. 



3 Some comparisons with other programs 

We conclude this report on the update of ZFITTER with few comparisons with other 
programs . 

For LEP 1, we restrict ourselves to the two small tables in Table 0. They represent 



an update of Tables 37 and 38 of |T^, which were produced with ZFITTER v. 5. 20. The 
complete Tables 37-40 (the latter two for another acollinearity cut) may be updated 
with use of our Tables |^ and ^. The agreement was not considered to be satisfactory 
in [|l^] but is excellent now. 



For a wider energy range, we performed a first comparison in Figure 3 of [|T3 



based on ALIBABA v.2 (1990) TOPAZO v.4.3 (1999) [|1|,|1|,^], and ZFITTER v.5.14 
(1998). We have repeated the comparison with the same version of ALIBABA, TOPAZO 
v.4.4 (1999) |T|,0, and ZFITTER v.6.11 (1999) in Figure 0. The ZFITTER numbers 



are produced with the default settings (if not otherwise stated), and the other two 
programs have also been run by ourselves. 

Concerning the TOPAZO v.4.4 ratios, we register a different behaviour (compared to 
v.4.3) for ^acc = 40° which is now much closer to the ALIBABA ratios. In ZFITTER we 
varied the treatment of some higher-order corrections via flags F0T2 and PAIRS with 
not too much effect. Furthermore, when we switched off the two- loop contributions in 
ALIBABA (with setting I0RDER=3), the agreement became much better. 

On the other hand, a cross check of the ZFITTER and TOPAZO programs applying 
s'-cuts comparable to the acollinearity cuts used for the figures above and below show 
a very high level of agreement between the two, at LEP 1 (< 0(3 ■ 10"'')), but also at 
LEP 2 energies at the order of less than a per mil ^ Initial-state pair production and 



sPor LEP 2 energies flag FINR was set to for the final-state corrections - recommended choice. 
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Tabic 4: Comparison of ZFITTER v. 6. 11 (first row) with ZFITTER v. 5. 20 (second row) for the muonic 
forward- backward asymmetry with angular acceptance cut (6'acc = 0°,20°,40°) and acollinearity cuts 
(^'acoi < 10°) and (^acoi < 25°). The initial-final state interference is switched off and only final-state 
radiation is corrected. 



exponentiation of higher-orders do not spoil this high level of agreement for the s'-cut. 
This may be seen in Figures [16 . 

Our s'-cut dependent ratios deviate from unity mostly in the regions where the 
radiative return is not prevented. The same is true for the ratios with acollinearity 
cut; since this cut is not as effective in preventing the radiative return as the s'-cut, the 
deviations survive at higher energies to some extent. This fact and the higher order 
corrections, which remained untouched by our study, seem to be the main sources of 
the remaining deviations between the different programs. 

In this context, it will be quite interesting to see the comparisons with s'-cut of 
ZFITTER and KORALZ |3^ and KK |3^ being extended to situations with acollinearity 
cut (see also 131,133). 
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Figure 10: Ratios of muon-pair production cross-sections and differences of forward-backward asymmetries 
predicted from ZFITTER v. 6. 11 and v. 5. 20 without and with acceptance cut and with three different 
acollinearity cuts: ^acoi < 10°, 25°, 90°, Emin = 1 GeV; programs differ by final-state radiation. 



4 Summary 

We derived analytical formulae for the photonic corrections with acollinearity cut and 
got substantial deviations from the coding in ZFITTER until version 5. The essentials 
of the changes have been described and numerical comparisons are performed in great 
detail. Fortunately, we may conclude that the numerical changes are not as big as 
one could expect. Although, certain differences to predictions of other codes remain 
untouched. They are pronounced when the radiative return to the Z resonance is 
kincmatically allowed. If one is interested to perform investigations in this kinematical 
regime further studies are needed. 
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with ^acol < 25° 


^acc 


M,-3 


- 1.8 


M, 


+ 1.8 




0° 


-0.28677 


-0.17066 


-0.00043 


0.11302 


0.15695 




-0.28647 


-0.17049 


-0.00043 


0.11293 


0.15682 


20° 


-0.27752 


-0.16514 


-0.00038 


0.10952 


0.15214 




-0.27722 


-0.16497 


-0.00037 


0.10944 


0.15204 


40° 


-0.24474 


-0.14562 


-0.00027 


0.09684 


0.13461 




-0.24445 


-0.14545 


-0.00026 


0.09678 


0.13454 



Table 5: Comparison of net corrections from ZFITTER v. 6. 11 (first row) with ZFITTER v.5.20 (second 
row) for the muonic forward-backward asymmetry with angular acceptance cut (^acc = 0°,20°,40°) and 
acollinearity cuts (^acoi < 10°) and (^acoi < 25°). The initial-final state interference is switched off. 
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Figure 11: Net ratios of muon-pair production cross-sections predicted from ZFITTER v. 6. 11 and v.5.20 
without and with acceptance cut and with three different acollinearity cuts: ^acoi < 10°, 25°, 90°; Emin = 1 
GeV; initial-final state interference not included. 
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Figure 12: Net differences of muonic forward-backward asymmetries predicted from ZFITTER v. 6. 11 and 
v.5.20 without and with acceptance cut and with three different acollinearity cuts: ^acoi < 10°, 25°, 90°; 
Emin = 1 GeV; initial-final state interference not Included. 
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Figure 13: Net ratios of muon-pair production cross-sections predicted from ZFITTER v.6.11 and v.5.20 

without and with acceptance cut and with three different acollinearity cuts: ^acoi < 10°, 25°, 90°; Emin = 1 
GeV; all the discussed corrections included. 
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0.31680 
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with ^acol < 10 


° 
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-0.14542 


-0.00011 
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0.13483 



Table 6: Comparison of net corrections from ZFITTER v.6.11 (first row) with ZFITTER v.5.20 (second 
row) for muon-pair production with angular acceptance cut (^acc = 0°,20°,40°) and acollinearity cuts 
(^acoi < 10°) and (^acoi < 25°). The initial-final state interference is switched on. 
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Figure 14: Net differences of muon-pair production asymmetries predicted from ZFITTER v.6.11 and v. 5. 20 
without and with acceptance cut and with three different acollinearity cuts: ^acoi < 10°, 25°, 90°; Emin = 1 
GeV; all the discussed corrections included. 
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with ^acol < 10° 
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Table 7: A comparison of predictions for muonic cross-sections and forward-backward asymmetries around 
the Z peak. First row is without initial-final state interference, second row with, third row the relative 
effect of that interference in per mil. 
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Figure 15: Muon-pair production cross-section ratios with 6'™^f = 10°, 25° and 6'acc = 40°; (a) TOPAZO 
v.4.3 and v. 4. 4 versus ZFITTER v. 6. 04/06 and v. 6. 11 (1999), (b) ALIBABA v.2 (1990) versus ZFITTER 
v.6.11. Flag setting: ISPP=0. 
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Figure 16: Comparison of predictions from ZFITTER v.6.11 and TOPAZO v. 4. 4 with s'-cut: (a) Muon- 
pair production cross-section ratios, (b) forward-backward asymmetry differences. Flag setting: ISPP=1, 
FINR=0; further: SIPP=S_PR |||. 
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A Radiative Return and Acollinearity Cut 

An acollinearity cut may act as a simple cut on invariant masses and thus it may prevent 
the radiative return of a/s' to the Z peak (and the development of the radiative tail) 
for measurements at higher y^. In Figure ^ it may be seen that a reasonable analogue 
of a cut value \J s"**" is the upper value of s' of region III, Re^^^i, defined in (|10]): 



Re 



with 



Rf) 



sm((?--/2) 



1 + sin • 
The relations are visualised in Table [A.l| . 
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0.4465 
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0.3333 
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Table A.l: The minimal cms energy, V s""", at which the radiative return to the Z peak is prevented by 
an acollinearity cut given as a function of this cut. 
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